Parametric amplification with a friction in heavy ion collisions 



Masamichi Ishihara * 
Department of Physics, Tohoku University 
Aoba-ku, Sendai 980-8578, Japan 
(February 1, 2008) 

We study the effects of the expansion of the system and the friction on the parametric amplifica- 
tion of mesonic fields in high energy heavy ion collisions within the linear a model . The equation of 
motion which is similar to Mathieu equation is derived to describe the time development of the fields 
in the last stage of a heavy ion collision after the freezeout time. The enhanced mode is extracted 
analytically by comparison with Mathieu equation and the equation of motion is solved numerically 
to examine whether soft modes will be enhanced or not. It is found that a strong peak appears 
around 267 MeV in the pion transverse momentum distribution in cases with weak friction and high 
maximum temperature. This enhancement may be extracted by taking the ratio between different 
, modes in the pion transverse momentum distribution. 
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I. INTRODUCTION 

t-H ! 

, The chiral symmetry restoration has been investigated theoretically and is expected to be realized in the coming 
experiments of high energy heavy ion collisions at RHIC and LHC. The chiral symmetry restoration leads to the 
Q\ [ subsequent rebreaking which has been studied in terms of the interesting phenomena called 'Disoriented Chiral 
f**«> ■ Condensates'(DCC). DCC is a (quasi) ground state and is characterized by finite pion condensates in the linear 
a model. Though DCC consists of zero mode in an ideal case, it consists of soft modes in a realistic case. In 
the conventional formation scenarios of DCC, it is assumed that the condensate rolls down to any direction with 
equal probability from the top of the hill of the mexican hat potential if the thermal equilibrium above the critical 
temperature is completely established. However, this is not always valid in heavy ion collisions at high energies because 
the initial condition is such that the condensate is settled on the minimum of the potential at zero temperature before 
the collision. We have shown that the condensate moves almost along the sigma axis when the effects of random 
forces are neglig ible H|). In such cases, the condensate moves to a (normal) or —a direction which is the opposite 
Cl^l side to the normal vacuum. DCC does not appear if the condensate stays at the normal vacuum, while if it stays at 
— cr, then it is just a DCC. 

As mentioned above, the condensate may oscillate along the sigma axis around the minimum of the effective 
potential at later time in high energy heavy ion collisions because the condensate moves along the sigma axis. The 
field which couples to the oscillating field may be amplified due to the so called 'parametric amplification' (^-||. It was 
pointed out that the enhancement of nonzero mode due to parametric amplification on the normal vacuum may occur 
when the condensate oscillates around the minimum of the potential at zero temperature . It was also shown that 
a squeezed state is formed by parametric amplification pj. However, it is not well-known whether the enhancement of 
nonzero mode occurs or not when the amplitude of the oscillation decreases because of the expansion and the friction. 

The equation of motion for the condensate with a friction term was derived by Biro and Greiner on the chirally 
symmetric vacuum where T > T c with T being the temperature of the system and T c the critical temperature, and 
by Rischke || on chirally broken vacuum where T < T c . The magnitude of the friction for the condensate on the 
chirally broken vacuum in the linear a model was given by Rischke. There is a finite friction for a field because of 
a — > 2tt decay. Then, it is better to include the friction in the calculation of parametric amplification for soft modes 
like DCC formation and decay studies |^-|lC|]. 

The aim of this paper is to investigate whether the enhancement of nonzero modes occurs or not in one dimensional 
scaling case when there is a friction. This paper is organized as follows. In sec.|n|, the equation of motion is derived 
and a constant friction is introduced to describe the dissipative nature effectively. After that, the equation of motion is 
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analyzed by comparing with Mathieu equation. In sec. Ill, The initial condition is given and a measure to extract the 
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amplification of the fields is introduced. The equation of motion is solved numerically to see whether the enhancement 
of soft modes occurs or not. Sec. IV is assigned for conclusions and discussions. 



II. EQUATION OF MOTION IN ONE DIMENSIONAL SCALING CASE 



We use the linear a model to describe the parametric amplification in the last stage (after freezeout) of high energy 
heavy ion collisions. The Lagrangian is 



V(cf>) = ±(ct> 2 -v 2 ) 2 -Ha>, 



where cj> = (a' , it). The minimum of the potential is determined by differentiating V (</)): 



da 1 



0. 



(1) 



(2) 



<r'=f„,n=0 



where f n is the pion decay constant. Since we are interested in the motion around the minimum of the potential, the 
new field a = a' — f v is introduced and the potential is rewritten. We have the relations: 



2 def . , , j 
"V = A(3/ ff - v ), 



«v = ML -v), 



H = \fMl-v 2 ), 



(3) 



where m CT is the mass of the sigma meson and m w is that of the pion. The equation of motion for a and 7Tj(z = 1, 2, 3) 
become 



Ua(x) + m\o(x) + A (<r 2 (x) + 3/ w cr(s) + tt 2 (x)) a{x) + A/^vr 2 = 0, 
□7Ti(ar) + mlmix) + A (a 2 (x) + 2f 7r a(x) + n 2 (x)) TT t (x) = 0, 



(4a) 
(4b) 



where □ = d 2 — d 2 — d 2 — d 2 . We solve the above equations of motion, especially for soft modes in the presence of the 
background oscillation of the condensation (zero mode) . Since the sum of hard modes generates a friction @J|,|ll[ in 
general, a friction term is introduced phcnomcnologically. The equations of motion become 

□ ct(x) + r] a (x)d T a(x) + m 2 a(x) + A (<r 2 (x) + 3f n a(x) + vr 2 (a;)) a{x) + Xf^n 2 = 0, (5a) 
□71,(0;) + r] 7T (x)d T n(x) + mln(?) + A (cr 2 (x) + 2f 7r a(x) + n 2 (x)) mix) = 0, (5b) 

where Tja-(x) and ^(x) are the frictions for a and tt fields, respectively. Convenient variables r and r\ arc defined by 

't + z 



def 



def 1 . 

V = 77 ln 



t - z 



for one dimensional scaling case. D'Alembertian is rewritten by using these variables: 

U=d 2 T + -d r -^d 2 -dl, 



(6) 



(7) 



where 9 2 = d 2 + d 2 . We assume that the frictions r) a and r? w depend on only r. First, we solve eq.(|a|) without the 
interaction term. The notation (r) which implies zero mode (the condensate) is introduced. It is the background 
field for nonzero modes. The equation becomes 



d A T +[vAr) + -)d T +m, 



<jW{t) =0, 



(8) 



where it is assumed that is a function of only r. Without loss of generality, one can substitute the following 
factorized form of <t( )(t) into eq.(H): 



a (o) (T) = / (o) (T)u (o) (r) _ 



(9) 



du^\r)/d,T term appears on the right hand side of eq.(||) and its coefficient vanishes if /^(r) satisfies the following 
equation: 
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' dr 



- + (^(r) + ;) / (0) (r)=0. 



(10) 



The equation for (r) becomes 

d 2 u^(r) 
dr 2 + 



J_ _ l dry CT (r) _ 1 
2r 2 2 dr 4 



(^(r) + ^ 2 +m4 u (°)(r)=0. 



(11) 



For a large r, 1/t and dr\ a jdT are negligible if r\ a reaches a constant. Then, eq.(|ll[) for large r is approximately 



^^ + {mj-^(oo)}«(?)(r)=0. 
The solutions of cqs.([To|) and (|l2| ) arc 

/(°)( r ) = Cr-^exp (-1 J T ds Va ( S )^ ~ Cr-^exp (-^(oo)r) , 



= tto cos 



77 2 (oo) 



(12) 

(13a) 
(13b) 



wher e C, uo and 9 are constants. Then, 0(A°) solution of g^\t) = /(°) (r)ui°' 1 (r) is given with the help of eqs.(13a) 
and dl3b| ). We can obtain the same expression for it fields. However, we now consider the case in which the motion 

of 7r field is negligible, namely 7r,( ^ ~ 0, because the condensate moves almost along the sigma axis. 

We solve eqs.(pa|) and (5b) with the interaction term in order to evaluate the effect of the motion of the condensate. 
Fields are decomposed to substract the background oscillation as follows: 



a (°)( x ) = aaT-^e-^A^T CO s 
n(x) = tt (0) (a?) + tt (1) (a:) ~ tt (1) (x), 



where ctq = Cuq which is a constant. With the help of eqs.(14al) and (14c), we have 



(>>) 2 + 2a^a^ + 



O 



JO) 



+ 2a( ( V 1 > + 0(A 2 ), 



it 2 = 0(A 2 ) 



(>>) + 3 + 3*^ ( CT d)) 2 + ( CT d)) 3 

a] +0(a(°)A 2 ) +0(A 3 ), 



O 



o 



TT 2 a 



0(A 2 ) x +0(A)) 



(14a) 
(14b) 
(14c) 

(15a) 
(15b) 

(15c) 
(15d) 



O ^[c^ -*] 2 ^ and higher contributions are negligible compared with 0(a^) because Co is taken to be small and 

t~ 1//2 cxp (— 77 CT (oo)r/2) is smal l for large r. The lowe st co ntribution of the interaction terms is obtained by neglecting 
0(A 2 ) and higher terms in eq. (15a), (15b), (15c) and (15d): 



3 2 + ( r h (T) + - ) d T + uoI{t) + 6\Ua^ (r) \ a W (r, k Tl k n ) = 



(16) 



where u^^t, kx, k^) is the Fourier transformation of cr*- 1 - 1 (t, x), kx is the transverse momentum, k v is the conjugate 
momentum of the 'rapidity' rj [eq.@] and w 2 (r, kx,k v ) = k T + k 2 /r 2 + m 2 which does not depend on k v for large 

t. erW (t, kx, k v ) — f(°^ {t)u^ (t, kx, k v ) is substituted into eq.(|l^) as eq.(^) is substituted into eq.(||) . Note that 
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/ (0) (r) is the same function as found in eq.( |13a| ) except for an overall constant. Then, the equation for at large 
r becomes 



d 2 Ua (r, kT,k v ) 



dr 2 

+6A/ w ct t" 



1 



l/2 e -i^(oo)r cos 



77^(00) 



1/2 



uW(r,*r,A,) = 0, 



where w^(fcr) = 00, fey, fe^). Applying the change of variable, 



2£ 



2 / \ \ V2 



to eq.(|17|), we obtain 



d 2 u^\^, k Tl k n 



de 

A <y {k T )=A 
QAO = - 



+ K(fc T ) - 2?,, (0 cos(20] u^K, fcr, M = 0, 

^(fc T )-^M/ 4 s 
™2 - 7/£(oc)/4 
12A/ X er 



r -l/2 e -^ CT (oo)7 



The field o^ 1 ) is obtained as follows: 

a«(T, fer,*,) = CrT-^e-^^uPiT, k T ,k v ), 
where C CT is a constant. In the same manner, we obtain the equation for tt^(t, kr, k v ): 



d^ij^K) + [MkT) 2qA0 cos{m ^ K) 



0. 



«*(0 = - 



^(fcr) -^(oo)/4 
777,2 - 772 (oo)/4 



4A/ W cr 



7--l/2 e -^<r(«))T 



where £r (for) = fej 



777-2. rp ne ^(1) j g a j go j)j- a ; ne( j ag f nows: 

7r« (t, fer, = C.T- 1 / 2 e -5^(-)- u ( r i)( Tj ferj 



(17) 



(18) 

(19a) 
(19b) 
(19c) 

(20) 

(21a) 
(21b) 
(21c) 

(22) 



where CL is also a con stant. The dependence of kr and k v is included in only A a and A v in eqs.( |19a| ) and (21a). 

Eqs.( [l9a|) and ( |21a| ) are similar to Mathieu equation which is the basic equation for the parametric amplification. 
Mathieu equation is 



d 2 F(Q 

de 



[A - 2qcos(2£)]F(£) = 0. 



(23) 



The coefficient in front of cosine is a constant in Mathieu equation, while it is a time dependent function, i.e. a 
function of £, in the present case. The stable and unstable regions of the solution of Mathieu equation are described 
in the q-A plane. The unstable regions at q ~ correspond to A — 1, 4, 9, • • •. The time dependence of q in the present 
case do not change enhanced modes corresponding to A = 1, 4, 9, ■ ■ •, while it changes the magnitude of amplification. 
In the linear a model, the friction coefficient of tt is almost zero at zero temperature, while that of a is not. The 
denominator of A a is not negative for r] a (00) < 277V. Then, the minimum of A a is not less than 4 , while that of A v 
is smaller than 1. Only ir has the amplification mode for A = 1. The amplified modes are determined by A = 1, 4, • ■ • 
and has 77 dependence. The 77 dependence is weak for r] 2 /{Am 2 ) << 1. 772/(47772) is about 0.1 for 77 CT = 2fm~ 1 and 
m<j = 600MeV. In this study, it is apparent that the 77 dependence of the amplified modes is weak since this condition 
is hold in the following calculations. Once the parameters of the linear a model and the magnitude of the friction 
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are chosen, the amplified modes are calculated easily. As an example, f n = 92.5 MeV, v = 87.4 MeV, A = 20 are 
chosen for numerical calculations in the next section. These parameters generate m a ~ 600 MeV, m n ~ 135MeV. The 
amplified modes for ri„ = 0.5fm _1 and -q^ = Ofm^ 1 are as follows. The amplified modes kr for 7r fields is about 267 
MeV for A„ = 1 , 583 MeV for A„ = 4 and so on. That for a field is about 669 MeV for A a = 9. A peak or peaks will 
appear in the kx distribution if the parametric amplification is strong enough. There is no mode which corresponds 
to A a — 1 and the mode for A a = 4 corresponds to zero mode (condensates). However, zero mode is the background 
field for the finite modes and a^\x) does not include zero mode. Therefore A a = 4 cannot be a candidate for the 
parametric amplification. 



III. NUMERICAL CALCULATION 



The remaining task is to solve eqs.(21a) and (19a) numerically with an adequate initial condition. The amplitude 
of the oscillator is given as an initial co nditio n. On the other hand, we choose some magnitude of the friction referring 
the results by some authors. In subsec.III A, the initial condition and the initial momen tum d istribution to solve the 
equation are chosen to mimic the situation in high energy heavy ion collisions. In subsec. III B| , we discuss the method 
to extract the enhancement due to parametric amplification in the transverse momentu m dis tribution. In subsec.III C , 
the equation of motion is solved numer ically with the initial condition given in subsec. III A and the characteristics of 
the enhancement formulated in subsec.III B is shown explicitly. 



A. Initial condition and magnitude of the friction 

In high energy heavy ion collisions, the system can be described roughly by scaling-hydrodynamics E^ ] in the 
cooling stage. The time dependence of the temperature in the cooling stage is 



T,, 



1/3 



T(r)=T m ^j 9{t } -t) for r > r m , (24) 

where T m is the maximum temperature of the system, r m is the time at which the temperature becomes maximum 
and Tf is the freezeout time at which the temperature becomes zero suddenly. If Tf is given, r/ is calculated as 
follows. 



r f = r m ^j . (25) 

Though the friction is a time-dependent function in general, it is expected that its magnitude reaches a constant 
because the temperature becomes (almost) zero after the freezeout. Then, we use the friction around the true minimum 
of the potential at zero temperature. The magnitude on the chirally symmetric vacuum at T ~ T c obtained by Biro 
and Greiner is 2.2 fm" 1 and the magnitudes on chirally broken vacuum at T = obtained by Rischke are about 
3fm _1 for a field and for 7r fields ||. The contribution from nonzero modes are summed to obtain the friction for 
zero mode in such calculations. However, the friction used in this paper comes from only hard modes. Its magnitude 
is smaller than that obtained by summing the contributions from all nonzero modes. For these reasons, the following 
values are used in the numerical calculations: 

X]„ = (0.25,0.5,1.0,2.0) far 1 , %=0. (26) 

There are other three parameters in the linear a model. We use the same values as in the previous section: = 92.5 
MeV, v = 87.4 MeV and A = 20. 

In the previous papers QJ^] , it was found that the initial condition is important for the chiral symmetry restoration 
in high energy heavy ion collisions. Then, it is conceivable that the initial condition for a heavy ion collision affects 
the initial condition for the parametric amplification. As already stated, the condensate moves along the sigma 
axis because it reflects the initial condition before the collision. The probable initial condition for the parametric 
amplification in heavy ion collisions is 



= (<7 ini , 0,0,0). (27) 

T = TJ 



The time derivatives of the fields are needed to calculate the time development. The condensate is almost near 
the minimum of the finite temperature effective potential in strong friction cases (of course, not too strong) and 
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the velocity of the condensate is small enough. On the contrary, the velocity is not small in weak friction cases. 
The amplitude of the condensate in nonzero velocity cases is larger than that in zero velocity case. The larger the 
amplitude of the condensate is, the stronger the amplification is. The amplification with the initial condition that 
the time derivatives of the condensate are zero is weakest compared with the amplification with the initial condition 
that the time derivatives are not zero. It will be concluded that the parametric amplification will occur in general if 
the amplification is found in zero velocity case. The velocity of the condensate at the initial time is thus taken to be 
zero. Then, the angle 6 in cq.(|l7|) is determined by da^ /dr\ T _ r ^ = : 

t 1 + ^(4/ \ ( 2 ^M xl/2 



9 = arctan '" v 11 U2 - roj. - r f . (28) 

The momentum distribution is proportional to the square of the fields. On the other hand, the momentum dis- 
tribution depends on only kx because of the scaling property. Therefore, the field which is the square root of the 
momentum distribution depends on only kx- Then, we take the initial value of a^'(rf] kx, k v ) and 77^(77; kx, k v ) as 

a^\Tf,kx,k v ) oc \] g{kr\m a ), 7T (1 ) (t/; k T , k v ) oc \J g(kx] m n ). (29) 

where g{kx', m) is the initial distribution function. The absolute magnitude which is not given in eq.(p9|) is unimportant 
in the following discussion because we will consider only the relative ratio. If the system is locally thermalized at the 
initial time, g{kx', m) is chosen as follows: 



CX P 1 — xj 



B. Extraction of the peaks due to parametric amplification 



One important problem is how to extract the peaks in the transverse momentum distribution generated by the 
parametric amplification in numerical calculations. The amplitude a and 7r become zero asymptotically for all mode 
because of the expansion and the dissipation. First, we introduce the following quantities to extract the magnitude 
of the amplification of the fields. 

1 a \ - u 1 u 1, \ dcf w ^ ( r ; fc T,fc„) rQ1 , 

R a (T]A a: Tf) = R a (T- l k T ,k ri ,T : f) = — -r , (31) 

i/ict (Tf,kx,k v ) 

where a = 7r, a and ip!)p = tt^\ tp^' = . The ratio R a oscillates as a function of r. The amplification of the 
fields is seen explicitly by removing the damping factor, exp[— r] a {oo)T /2]t^ 1 / 2 from R a . It is performed by using the 
following quantity: 

r a (r;A a ,r f ) d ^ ui 1 \r,k T )/u^(r f ,kx). (32) 

It is useful to consider the envelop of R a and r a for our purpose. The envelope of r a is denoted by r e a . However, 
R a is not an adequate quantity to extract the amplification and the envelop of R a too because of the damping. 
Fortunately, the damping factor is momentum independent. The reduction of the amplitude by the expansion and 
the dissipation are cancelled by taking the ratio between different modes. We use the r e a because the envelop is easily 
calculated for r a . 

P a (A a ,A a , r ) d ^ lim r ; f iT; f aiTf \ , (33) 

where A^^ and A a r are the reference to other modes, A v and A a respectively Since the ratio of the momentum dis- 
tribution is proportional to [-P(^4, ^4 r )] 2 , it links directly the amplification and the transverse momentum distribution. 
The effects of the initial distribution are taken into account by multiplying the square root of the initial function as 
follows: 



G 



P«™(A,A r )^P(A,A r )Jj^, g(A)=g(k T ), (34) 

where g(kx) is the initial distribution function introduced in eg. (|30|) . The transverse momentum distribution can be 
calculated by eq.(|34|) for various initial distributions. 

The time derivative terms are not explicitly included in the quantities introduced in this subsection because the 
effect of these terms are small. These effects in terms of above quantities are evaluated in the appendix 
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C. Numerical Results 



What we would like to know are 1) whether parametric amplification occurs or not, and 2) whether its signal can 
be observed experimentally or not. The answer to the question 1) can be obtained by studying r a . Hereafter, we 
consider only the 7r fields because a field will not be amplified strongly. The reasons that the sigma field is not 
amplified strongly are that 1) the smallest amplified mode corresponds to A a = 9 , and 2) the effect of the friction 
for a |eq.(pO|)l does not vanish, while that for it vanishes [eq.(^2])] a t zero temperature. 

First, we show the time development of and r„ for A n = 1, <7i n ; = — 30MeV, r/ a — 0.25fm _1 and Tf = lOOMeV. 



The evolution of is shown in FigJld and that of r n is shown in Fig.fb, The amplitude of R^ increases temporarily 



because of the parametric amplification and then decreases because of the expansion and the dissipation. Removing 



the effect of the expansion and the dissipation, we can see the amplification obviously shown in Fig. lb. The amplitude 



of r w becomes a constant since the factor q w (£) [see eq.(21c)] becomes zero. 

The asymptotic value of the amplitude, r^(r = oo, A n ), is shown in Fig.|| for various T m as a function of A^ 
for (Tini = -30MeV, rj a = 0.25fm _1 and T f = lOOMeV. There is a peak near A v ~ 1 caused by the parametric 
amplification. This strong amplification may give the possibility of the observation even when the initial distribution 
does not have enough amplitude at A ^ 1. It seems that the amplification saturates as T m increases. This is 
understood as follows. The increase of the amplitude by the parametric amplification depends on the initial amplitude 
of the condensates and the time interval (s = r —Tf). Since the amplitude at the initial time is fixed in this calculation, 
the time interval is investigated here. Amplification described in Mathieu equation is characterized on the q-A plane, 
where q and A appear in eq. (p3|) . The field is amplified if the set of the parameters q and A is in a unstable region 
on the q-A plane. The set (q, A) at the initial time is shown by dot in Fig||. The coefficient g(£) defined by eq.(21c) 



is a time dependent function and it moves to zero as shown in Fig|l, while it is a constant in Mathieu equation. The 
point on q-A plane goes out of the unstable region as time goes by. This implies that the small est q(£) for each modes 
exists for the field to be amplified. Then, we consider a certain time tq at which q v (^) in eq.( |21c| ) becomes a certain 
value, Q. We consider also Q/(t/) which is when r = Tf. With the help of eq.(|21c[), we obtain 



cxp (-vAoo)s) = I 1 + — ] , ( V>) 



Q}(Tf) V Tf 

where s = tq — Tf. The left hand side is a function of only s, while the right hand side depends on both s and Tf. 
That T m goes to infinity corresponds to that rt goes to infinity. Then, Tf — > oo lim it can be taken with a fixed s. 



Since Qf{jf) converges at a certain value Qf which is obtained from eqs.(14b), (21b) and (|28|), the right hand side of 



eq.(|3j|) becomes Q /Qf. Then, s is finite as T m becomes infinity, which is obtained by solving eq.(|E 

2 



Vt(oo) 



ln(Q/Q f ). (36) 



This implies that the time interval for the amplification is still a finite constant even when T m is infinite. We conclude 
that the amplification saturates as T rn becomes infinite because the initial amplitude is finite and the time interval 
converges. 

The parameters, 7/0- (oo), (Ti n i and Tf, are varied to see the sensitivity of the amplification in the following calculations. 

The ratio of the amplifications, r%{r — 00, A n ), is shown in FigJ^ for various r/ a (oo) as a function of A n . The 
magnitude of the amplification decreases as the magnitude of the friction increases. It is reasonable because the 
oscillation of the condensate are reduced rapi dly wh en t he friction is strong. It is apparent that the amplitude of the 



oscillation decreases since 2g Q (£) in eqs.(21a) and (19a) vanishes quickly in the strong friction cases. It was found 



that the parametric amplification cannot occur for frictions which is larger than 2fuT 

The ratio of amplifications , r£(r = 00, A n ), is shown in FigJ| for various a ml as a function of A„. The dependence 
on <7i n i is important for the parametric amplification in the presence of the expansion and the friction because the 
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amplitude of the oscillator (condensate) decreases as time goes by. However, we can find a peak near ~ 1, at least, 
for Cini < — 10 MeV. Then, it may be possible to find this peak in the transverse momentum distribution. 

The ratio of amplification , r^(r = oo, A v ), is shown in Fig.^ for various Tj as a function of A^. It is obvious that 
Tf dependence is weak. 

Hereafter, the A„ distribution, namely the transverse momentum distribution, is displayed ta king into account the 



initial transverse momentum distribution. The thermal distribution [eq.(30)] is assumed and eq. (p4|) is used. We take 
A r = 0.5 which is the reference in p now defined in eq.(|34|). 

The ratio between different modes including the effect of the initial thermal distribution, P^ cw (A 7r , A w>r — 0.5), 
for various T m is shown in Fig.0. The parameters are a ini = — 20McV, Tf = lOOMcV and r? — O.Sfm -1 . The 
ratio between the momentum distributions corresponds to |P" ow | 2 . Then, it is found from Fig.0 that the ratio at 
An ~ 1 is 4 times larger than that at A^ = 0.5 for T m — 300MeV. The increase in small A n and the decrease in 
large A^ reflect the initial thermal distribution. Then, P ncw (A,A r ) for large A limit becomes zero when the initial 
distribution is thermal. Similar calculations can be performed to obtain P" cw for various parameters. For example, 
P° ew (A w - 1, A^ r = 0.5) - 12.7 for <j ini = -30.0MeV, T m = 300MeV, T f = lOOMeV and r] a = 0.25fea _1 . 

There are many A's at which a field is amplified in Mathieu equation. The amplification around A^ = 4 is also 
investigated, but strong amplification is not found. This comes from the property of Mathieu equation and the small 
initial amplitude in the initial thermal distribution compared with that at A^ ~ 1. Then, it is difficult to observe the 
effects of amplification except for A w ~ 1. 

IV. CONCLUSIONS AND DISCUSSIONS 

We investigate the effects of the expansion and the friction on the parametric amplification in one dimensional 
scaling case within the linear a model. The equation of motion at zero temperature is derived in the case where the 
sigma condensate oscillates around the minimum of the potential. From this study, we found that the amplitude of the 
oscillator (zero mode) and the magnitude of the friction are essentially important for the parametric amplification in 
one dimensional scaling case. In a weak friction case, the strong enhancement will occur at A v ~ 1. (fc-r ~ 267MeV) 
in the rjion transverse momentum distribution and it may be possible to be observed by the normalized distribution, 
pnowj _ Then, we have a chance to catch the signal of the temporal restoration of chiral symmetry [jj]. Contrary 
to this, it may be difficult for the field to be amplified by the parametric amplification in a strong friction case. On 
the other hand, high maximum temperature and/or sufficiently strong friction are needed for chiral symmetry to be 
restored for a long time as we found before pi . 

High maximum temperatures (above 100 MeV) was also assumed in this study, which settles the condensate at an 
unstable point on the potential. However, the mechanism of parametric amplification does not always require high 
temperature. The partial chiral symmetry restoration may be also observed if the amplitude of the condensate is large 
enough. 

It is possible to devide the modes of the field into three parts, the condensate, soft modes and hard modes. (See 
Fig.||) The hard modes can be regarded as the cause of the friction for the condensate and soft modes. There are the 
energy flows from the condensate and soft modes to hard modes. In the same way, there is the energy flow from the 
condensate to some soft modes. The decreasing of the zero mode amplitude comes from the energy dissipation due 
to the friction and the effect of the expansion. However, the energy flow from the condensate to the soft modes is not 
taken into account in the present study. Back reaction is discarded too. They will be taken into account in the future 
study. 
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APPENDIX A: THE EFFECTS OF THE TIME DERIVATIVE TERMS 

In this appendix, the effect of the time derivative terms of the fields is checked for the quantities introduced in 



[II B. First, we consider the effect on R a . The ratio of the number-density like quantities is evaluated in order to 



make it clear that the effect is small in the present case: 
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R^(r;A a ,T f ) 



dip^/dr 



^\ Tf ;A a ) 



~ -2 
+ W Q 



dip^/dr 



T=T; 



(Al) 



Note that Cj a does not have kf /t 2 term. It is ignored because r is large enough. The time derivative of ipa is 



#1^ 



exp (-r) a r/2) du a r\ a exp {-r) a T/2) 



1 exp {-r] a T/2) 



dr t 1 ' 2 dr 2 rV2 " 2 r 3 / 2 

where ry a is the friction for ipa . At the time (r e ) at which du a /dT = 0, one has 

7/ Q 1 \ exp(-?7 Q T e /2) 



r/r 



2r e 



(r e ) 



1/2 



(A2) 



(A3) 



" (2) 

The same expression is obtained at the initial time for the intial condition, (du a / 'dr) T _ = 0. R a at r e is evaluated 
since the amplification is evaluated at the envelop of r a . 



exp [-r] a (r e - tj)] 

T e/Tf 



l + {2u a )- 2 (y a +T- 1 ) 2 



2 def exp -Tfo (r e - Tf 2 
^ = ; D a {T e ,Tf,r]) r a . (A4) 



>J T f 



The effect of the time derivative terms is included in D a (T e ,Tf 1 rj). Apparently, D a (oo, Tf, rj a ) < D a (r e ,Tf ,ri a ) < 
D a {Tf,Tf,ri a ) = 1. Then D^{oo^ tj, 0) for tt^ and the minimum of D CT (oo, tj, rj a ) with respect to rj a for cr^ 1 ^ are 
evaluated with a fixed tj : 



A r (oo,r /> 0) = [1 + 1/ (4r|^)] 1 

D a (<X,Tf, fy)^ = 



16w2 T/ 2 + 1-1 



16d;2T y 2 + 1 + 1 



(5a) 
(5b) 



_D Q at fcr* = is the minimum on the kx axis. Tf is some 10 fm, m n is 135 MeV and m CT = 600 MeV in the present 
case. In such a case, and D a at the minimum are about 0.9947 and is about 0.9837 respectively for Tf — 10 fm 
and fc<r = The maximum difference between D a with time derivative terms and without is less than 2%. These are 
small enough. 

The effect of the time derivative terms for P a (A a , A a ^ r ) is evaluated in the same way. P a {A ai A a ^ r ) is the ratio 

~ (2) 

between different modes at the same time, while Ra is the ratio between different times at the same mode. It is 
shown that time derivative terms do not affect on P n , while it is shown that time derivative terms affect slightly on 
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FIG. la. Time development of the ratio with A n = 1, a ini = -30MeV, T m = 300MeV, T f = lOOMeV and 
?? CT (oo) = 0.251m- 1 . 



FIG. lb. Time development of the ratio u n with A„ = 1, cr ini = -30MeV, T m = 300MeV, Tf = lOOMeV and 
77 CT (oo) = 0.25ftn _1 . 



FIG. 2. The amplifications of the field amplitude in the broad range T m . The initial condition is Ui n i = — 30MeV, 
t? ct = 0.25ml- 1 and T f = lOOMeV. 

FIG. 3. Sketch of the stable and unstable regions around A = 1 in the q-A plane. The shaded region is corresponding to 
the unstable region. The initial values in the q-A plane is indicated by black balls. The set (q, A) moves on the thick line to 
the q = direction with a constant A. 

FIG. 4. The amplifications of the field amplitude for various magnitudes of friction. The initial condition is T m = 300MeV, 
o-ini = -30MeV and Tj = lOOMeV. 



FIG. 5. The amplifications of the field amplitude for various ow. The initial condition is T m — 300MeV, r\ a = 0.25fm 
and Tf = lOOMeV. 



FIG. 6. The amplifications of the field amplitude for various Tf. The initial condition is T m = 300MeV, r\ a = 0.5fm 1 and 
o-ini = -30MeV. 



FIG. 7. The ratio of amplification of the amplitude for various T m with Ui n i = — 20MeV, Tf = lOOMeV and r\ = 0.5fm . 

FIG. 8. Sketch of mode division, ij and r)' are frictions. The arrow with the word 'energy' implies the energy flow. The 
effect of T) is ignored in the present study. 
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